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Abstract 

The coefficient functions for heavy-flavour production in deeply inelas- 
tic electron-hadron scattering have been calculated previously. Analytic ex- 
pressions are impossible to publish due to their length. Therefore we have 
tabulated them as two-dimensional arrays as is often done for the scale- 
dependent parton densities. Using this computer program we present event 
rates for charm production at HERA in bins of x and Q^. These rates are 
insensitive to variations in the factorization and renormalization scale /i. 



In the past few years, calculations of 0{as) QCD corrections to heavy- 
flavour production have made great progress (for a recent review see [|l|). 
Calculations have been completed for hadron-hadron collisions 0|,[3> P^^^- 
toproduction 0,0, electroproduction 0, 0|) and photon-photon collisions 
(real as well as virtual photons) [§ , • The calculation of higher order cor- 



rections to these processes are important for the the top quark search [|I0 
and the determination of the gluon distribution function |Tl[] , [^ which can 
be measured in open charm production. 

The expressions are only available in large computer programs for the 
radiative corrections as their complexity prohibits publishing them in an 
analytic form. This complexity is due to the non-zero mass m of the heavy 
quark. If m were zero the final formulae could be constructed and published 
as is shown e.g., for the 0{a1) corrections to the coefficient functions in 
deeply inelastic lepton-hadron scattering [|T^] and the Drell-Yan process |I^]. 
In some special cases the use of lengthy computer programs can be avoided 
by either making approximations [|15| or by making algebraic fits to the exact 
coefficient functions of heavy-flavour production in hadron-hadron collisions 
[|] and in photon-hadron scattering H. Such fits are quite accurate and 



enabled the authors in [^ to present complete tables of cross sections for 
charm, bottom and top production in many hadron-hadron reactions. 

These algebraic fits mentioned above could be made because the coeffi- 
cient functions calculated for heavy-fiavour production in hadron-hadron and 
photon-hadron collisions only depend on one scale independent variable. This 
is in contrast with the coefficient functions of deeply inelastic heavy-flavour 
production which depend on two scale independent variables as the photon 
is virtual. Moreover it turns out that in the latter process the coefficient 
functions show a much more complicated behaviour than in the former ones 
so that an algebraic fit as presented in is very difficult to achieve. There- 
fore we have constructed a set of tables reproducing the coefficient functions 
for each partonic subprocess in deeply inelastic electroproduction of heavy 
flavours. These tables are presented in the form of a two dimensional array in 
a computer program, analogous to the way various groups |T^, |]IB[ present 



the scale dependent parton densities in hadrons, which also depend on two 
variables. The computer program is available. Q Our two-dimensional tables 
are used to predict the 0{as) corrected rates for inclusive charm production 

"'^ Requests should be sent to smith@elsebeth.physics.sunysb.edu. 



at HERA by adding integrations over bins in x and Q^. It will be interesting 
to compare our results with the forthcoming data from the ZEUS and HI 
collaborations. Notice that in deeply inelastic electroproduction we only have 
to contend with the parton densities in the proton. In a later paper we will 
address the complications associated with predicting rates for heavy-flavour 
photoproduction (small Q"^) which also involve the parton densities in the 
photon. 

We begin by listing several important formulae. Omitting charged-current 
interactions, heavy-flavour production in deeply inelastic electron-proton scat- 
tering proceeds via the reaction 

e^ih) + P{p) ^ e-ik) + Q{pi){Q{p2)) +X. (1) 

Here X stands for any hadronic final state allowed by quantum number 
conservation. We sum over these states so that the process is inclusive with 
respect to the outgoing hadrons. When the virtuality — g^ {q = h — h) of 
the exchanged photon is not too large (— g^ <^ M"^) the reaction in (1) is 
dominated by the one-photon exchange mechanism and we can neglect any 
weak interaction effects. If we also integrate over the heavy (anti) quark Q 
(Q) in the final state the deeply inelastic electroproduction cross section can 
be written as 

"^'"^ ^''"' 5[{1 + (1 - yr}F,{x, Q^ m') - y'F,{x, Q^ m')] , (2) 



dxdy Q^ 

where 5* denotes the square of the cm. energy of the electron-proton system. 
The variables x and y are defined as 

^ = ir— (0<a;<l) , y = ^ (0<i/<l), (3) 

Ip ■ q p ■ h 

with 

-q' = Q' = xyS. (4) 

The deeply inelastic heavy-flavour structure functions appearing in the cross 
section (2) are given by ^2(0:, Q^,m^) and Fj,(x, Q^,m^) (longitudinal). The 
structure functions are given by the formula (see (6.5) of |^) 



Air'^m'^ Jx z 



2/^,2 PZ^^^ ^^ ^ nr> , , , , 1,2 
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Tim? 



r--^\2 r.X 2w (1) (1), /^^ 



'^ "- ^1^ ^(l)l^ /^^ 



+eL/4,/i^)(<;+<]ln£^)]], (5) 

where k = 2, L and the upper boundary on the integration is given by 2;max = 
Q'^/iQ^ + 4m^). Further fi^x^iJ?) , (i = g,q,q) denote the parton densities 
in the proton and /i stands for the mass factorization scale, which has been 
put equal to the renormalization scale. The coefficient functions, represented 
by c^lir], , dliv, . (^ = ^ . ? . ^ ; ^ = 0, 1) and by 4i(?7, , 4'k^ 0. (^ = 
q ,q;l = 0,1) are calculated in ^ and they are represented in the MS scheme. 
Furthermore they depend on the scaling variables t] and ^ defined by 

Am'' m^ 

where s is the square of the cm. energy of the virtual photon-parton sub- 
process which implies that in (5) z = Q'^/{Q^ + s). In this equation we made 
a distinction between the coefficient functions with respect to their origin. 
The coefficient functions indicated by Cki{'r],^),Ck{{'ri,^) originate from the 
partonic subprocesses where the virtual photon is coupled to the heavy quark 
whereas the quantities dki{T],$,) ,dj^{{7],^) come from the subprocess where 
the virtual photon interacts with the light quark. Hence the former are mul- 
tiplied by the charge squared of the heavy quark e%, whereas the latter are 
multiplied by the charge squared of the light quark e| respectively (both in 
units of e). Although terms proportional to ch^l appear in the inclusive 
photon-parton differential distributions they integrate to zero in the total 
partonic cross section, so we have not included them in (6). Furthermore we 
have isolated the factorization scale dependent term ln(/i^/r?7,^). The func- 
tions multiplied by this term, which are indicated by a bar, are called mass 
factorization parts. Notice that in the subsequent equations we discuss the 
transverse coefficient functions indicated by the subscript T instead of the 
ones indicated by the subscript 2. The relation between them is given by 

cg(^,0 = 4%0 + c%{v,0 and 4k^,0 = d%iv,0 + d%iv,0- where 
the same definition holds for the coefficient functions indicated by a bar. 
In the limit C, —>■ (see (6)) where the virtual photon becomes on-shell the 



above coefficient functions tend to their analogues obtained for photopro- 
duction which can be found in @],0]. For that purpose we had to modify 
the original expressions for the functions (i^^^(?7,,^). Originally c4,g('7'0 did 
not exist because no mass factorization was needed and in the limit ^ ^ 
d)p Jr],^) diverged logarithmically. This is due to an additional coUinear di- 
vergence which appears when the virtual photon coupled to the light quark 
goes on- mass-shell. Furthermore in the same limit <ii„(?7,0 did not vanish. 
Therefore, to use these functions in the region ^ ~ one has to perform 
an additional mass factorization to remove the coUinear divergence which is 
due to the on-shell photon. This is achieved by subtracting a term which 
is multiplied by a scale invariant function called R{^) = exp(— 20^) ( see 
(5.10) in ) where the subtraction is imposed if Q^ < Qfnin = 1-5 GeV^. 
This leads to the appearance of the function (i^^ ( rf}^ ^ = ) which would 
not have been present when the photon is treated to be highly virtual. This 
implies that the function dj^ will be proportional to -R(0 ^ind it vanishes 
when ,^ ^ oo. Notice that our choice of i? in was not scale independent 
so that the plot 11. b in that paper still contains a scale dependence. This 
however has no consequence for any of the numerical results. The above 
procedure implies that in the on-mass-shell limit the function ci^^ (77, 0) = 

and djpjj], Q),dr].J^'r], 0) become equal to the on-mass-shell photon coefficient 
functions in H, '|| (see (2.11) and (2.15) in §). 

The coefficient functions for the Born reaction (virtual photon-gluon fu- 
sion, see [|19|) are given by 



2 ' (1 + r; + ^/4)3 I- > "• '" (i + ,,)i/2 _ ,,1/2 






(7) 

7] \l/2 






eV , , o 1, (1 + 7^)1/2 + ^1/2 



where Tf = 1/2 in SU{N). 

No such simple analytic expressions can be given for the next to leading 
order coefficient functions. Therefore we present them in tables constructed 



as follows. First we ran the programs in and computed the coefficient 
functions for a grid of values of t] and ^ as defined in (6). This we have 
only done when they are represented in the MS scheme. We then divided 
these functions by the appropriate colour factor and then subtracted the 
asymptotic and threshold dependences, for which analytic expressions are 
available in the literature and will be presented below. Finally we wrote 
subroutines to set up two dimensional arrays. The interpolation is done in a 



bilinear fashion 21 



Starting with the virtual-gluon subprocess we define a new function with 
the threshold and asymptotic behavior removed, namely 



and 



h%,,{v, = (CATf)-' 4L{v, - PE.AV, ■ (10) 



Here we have split the coefficient functions c^ according to their colour parts 
indicated by the subscripts A and F. The colour factors are given by C^Tj 
and CpTf respectively, where for SU{N) ,Ca = N and Cp = {N^ - l)/2N. 
Further we have defined 

The mass factorization parts Ck,g(r],^) can be parameterized in a similar 

way by 

h'kliv,0 = {CATfr'cP^iv,0-PG,{v,0-pEkAv,0- (12) 

The functions E/^^Cy E^c with k = T ,L and C = A ,F describe the 
threshold behaviour as 77 — *> (or as s — i> Arn^) and are derived from (5.7)- 
(5.9) of IP] . The asymptotic behavior which holds in the region rj —>■ 00 
(or as {s —>■ 00) is given by the functions Gk,Gk with k = T ,L. The latter 
are obtained from [^ (see their Appendix A). The functions describing the 
threshold region have the following form 



^Notice that an extra factor of two should be muhiphed to (5.9) of 
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1 
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vr 
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6vr(l + e/4)= 



/3^ -ln(4/?2 



/3 -ln(4/5^ 



47r 1 + ^/4 
The functions describing the asymptotic region are given by 
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+ (l-^-i- ^ 



e 61 + e/4 
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where the functions J(^) and /(,^) are defined by (see Appendix A in ||2 



AO 



(4 + ^)1/2 + ^1/2 

(^(4 + 0)^/2 "^ ^ (4 + 0^/'-^/' 



In 



(23) 



TT 



' 1, 2({^+if"+e'' 



In 



(e(4 + 0)^/^L 6 2 V (4 + 01/2-^1/2 
..nH^ii#;^).2L,(<li#4z|^)l, (24) 



/(O 



2(4 + 0^/^ ^ '^ 2(4 + 0^/^ 

where Li2(a;) is the dilogarithmic function defined as 

[■X (U 

Li2(x) = - / -ln(l-t). (25) 

JO t 

We proceed in an analogous way for the coefficient functions corresponding 
to the virtual photon-light quark subprocesses. When the photon is coupled 
to the heavy flavour they are parameterized as 

h%,,{n.i) = {CFTf)-'c^^l{v,0-P'G,{v,0, (26) 

and 

h^Sav, = {CpTfY' c^j^liv, - P'G.iv, ■ (27) 

When the photon is coupled to the light quark we get 

hiUv,0 = iCFTf)-'d'j}l{v,0, (28) 

and 

h^lLiv,0 = iCFTfr'dil{v,0- (29) 

The subscripts H and L in the above expressions indicate that when they 
are inserted in (5) they have to be multiplied by the charge factors ejj and 
e| respectively. 

Notice that for the reactions discussed above only the A-part of the coef- 
ficient functions show an enhancement in both the threshold region and the 
asymptotic region. For the F-part we only observe large corrections in the 
threshold region except for the process given by the expressions in (28) and 
(29). To illustrate the quality of the fits, we present the plots of the coeffi- 
cient functions which constitute the bulk of the 0(0;^) radiative correction. 
They are given by hA^j.^^ (9) and hp^rpg (10), which are shown in fig.l and 
fig.2 respectively . They are plotted versus 77 for several different values of ^. 



(One can compare them with the plots in |^ for the corresponding functions 
without the subtraction of the threshold and asymptotic behavior.) In both 
instances, we observe that the h functions tend to zero as 77 — ;> and as 
7] -^ 00. Also, as ^ — i> 00, the h functions become zero across the entire 
77-region. We also see the complicated behaviour of the functions in the in- 
termediate region which illustrates why an algebraic parameterization is very 
difficult to accomplish. 

Using the fits, we now present single-particle inclusive event rates for 
inclusive c production given an integrated luminosity of 300 nb~^ at yS = 
298 GeV. Notice that we consider charm production only and do not sum 
over the charm and anti-charm cross sections. We take rric = 1.6 GeV/c^ and 
vary the mass factorization scale /x = y/Q'^ + m'^ up and down by a factor 
of two. Furthermore we have chosen the CTEQ2M parton densities |jT^ and 
the two-loop running coupling constant with Aq^j^ = 213 MeV. Our results 
are listed in tables 1 and 2 where we computed the cross section in (2) in 
bins of X and Q^ (these are the bins used by the ZEUS collaboration in their 
1993 data for F2). Table 1 contains the event numbers for the lower Q^ bins 
and Table 2 the corresponding numbers for the higher Q^ bins. 

We find the events concentrated at low Q^ and x, with an approximate 
five percent uncertainty at small x coming from the variation in the scale /i. 
In the intermediate Q^ and x region, we find the number of events stable and 
significant. As Q^ increases, the number of events drops rapidly. However, 
here we expect weak interaction effects to reduce the applicability of the 
one-photon exchange approximation. To conclude, we see that the number 
of charm events is large and relatively insusceptible to variations in the scale. 
From these results and the good luminosity at HERA, the extraction of the 
gluon density with significantly reduced uncertainty should be possible. 
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Figure Captions 

Fig.l. hArp,giv^O (^Q- (9)) versus r^ for ^ = 10~^ (solid line), ^ = 1 (dotted 
line), C, = 3.16 (short-dashed line), ^ = 31.6 (long-dashed line) and 
C, = 316 (dot-dashed line). 

Fig. 2. hp^j.g(r],^) (eq. (10)) versus r] hr ^ = 10"^ (sohd line), ^ = I (dotted 
line), C, = 3.16 (short-dashed line), C, = 31.6 (long-dashed line) and 
C, = 316 (dot-dashed line). 
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